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Abstract 

We analyze neutrino oscillations in a thought experiment in which neutrinos are produced 
by electrons on target nuclei. The neutrinos are detected through charged lepton production 
in their collision with nuclei in detector. Both the target and the detector are assumed to 
be crystals. The neutrinos are described by propagators. We find that different neutrino 
mass eigenstates have equal energies. We reproduce the standard phase of oscillations and 
demonstrate that at large distance from the production point oscillations disappear. 
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1 Introduction 



Neutrino properties are of great interest for particle physics, astrophysics and cosmology. 
One of the main sources of information about neutrinos is an investigation of neutrino os- 
cillations. Although a great many of papers on neutrino oscillations were published, there is 
still no common point of view on this phenomenon. One way to study neutrino oscillations 
is to consider neutrinos in the framework of the standard plane-wave description. Such an 
approach neglects effects which concern production and detection of neutrinos and demands 
choosing between two scenarios: 

1) equal momentum scenario (see papers by Gribov and Pontecorvo and by Fritzsch and 
Minkovski 0), 

2) equal energy scenario (see papers by Lipkin 13 and Stodolsky 4 a as well as by Kobzarev et 
al. [3] and Grimus and Stockinger 6 ). In ref.[Sj neutrino was created in /3-decay of a neutron 
localized at point P and detected through its interaction with an electron localized at point 
D. In refs.JS] and |Zj the authors considered neutrino produced by an electron beam on the 
target nucleus A and detected through its interaction with the nucleus B of the detector. 
The whole process looked as (see Fig. 1) 

e + A + B -> l + C + D, (1) 

where I is a lepton (e, \i or r), while C and D are recoil nuclei. In ref. j7j the nuclei A and B 
were supposed to be unconfined in a gaseous target/detector and described by wave packets; 
the electron wave function was also assumed to be a wave packet. 

In this article we investigate a thought experiment similar to that considered in ref. [7| but 
with nuclei A and B bound in crystals. We use the rigorous quantum field theory approach 
(Feynman diagram with neutrino propagators) to achieve the following goals: 

1) to show that in the case under consideration equal energy scenario takes place, 

2) to reproduce the standard form of the oscillation phase, 

3) to integrate over the phase space of the final particles in order to obtain the probability 
of the process, 

4) to demonstrate that oscillations disappear at large distances from the production point, 

5) to investigate corrections due to non-zero temperature and to show that they are not 
essential for the range of temperatures at which crystals may exist. 

One of the reasons to write this paper is that in the year 2004 in the most authoritative 
particle physics review of Particle Data Group the so-called " equal momentum scenario" was 
chosen to describe neutrino oscillations B.. We think, following Vysotsky jjjj, that though this 
approach gives the standard result, it is not self-consistent and thus misleading: neutrinos 
produced at point A (by electron) would not have a definite (electronic) flavor, but their 
flavor would oscillate with time at point A. 

In the case of solid state detector with stationary nuclei the neutrino energy is determined 
by energy conservation for the interaction in the detector and does not depend on neutrino 
masses. This conforms with the point of view of Lipkin 3 . Contrary to that the momentum 
cannot have a certain value because of the uncertainty relation for spatially localized nuclei. 

The paper is organized as follows. In section 2 we derive the amplitude for the process (^Q) 
and obtain the oscillation phase. In section 3 we integrate modulus of the amplitude squared 
over the momenta of the final nuclei and over the energy of the final lepton. In sections 2 
and 3 the temperature is supposed to be zero. In section 4 the case of non-zero temperature 
is considered. Our main results are summarized in section 5. Appendix contains derivation 
of factors which suppress oscillations at large distances. 
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Figure 1: Little donkey diagram from ref. [7j 

2 Nuclei bound in crystals at zero temperature 

In this section we follow the lines (and the notations) of ref. [7]. The difference is in the initial 
wave functions of the target and detector nuclei: they are stationary in this article while in 
ref.jTj they were represented by wave packets. 

Consider neutrino production by electron e on a target nucleus A of mass Ma- Then the 
neutrino collides with nucleus B of mass Mb in a detector and produces charged lepton I 
(see Fig. 1) . 

The electron neutrino v e produced on nucleus A is a superposition of three neutrino mass 
eigenstates: v e = Y2j Uej^j, where Vj is a state with mass rrij, U is a unitary mixing matrix, 
the first and second indices of which denote flavor and mass eigenstates, respectively. The 
detection of the neutrino by means of its interaction with B results in the projection of three 
propagating neutrino states onto the final flavour state v\ = ^ ■ Uiji/j. 

The amplitude of the process (0) is given by the following equation: 

A e — yl = ^ ^ UgjZJljAj , 

j 

where Aj is the amplitude for a given neutrino state of mass rrij 

Aj = J d 4 X 2 d 4 X 1 ^B(x2)^(x 2 )^Ux2)G j (x 2 - X^Aix^eix^c^l) ■ (3) 

Here ^f a are wave functions of particles a; Gj is the Green function of j-th mass eigenstate 
of neutrino, and x& = (ifc,x/%) are 4-dimensional coordinates. 

There exists a vast literature in which not only incoming but also outgoing particles are 
described by wave packets (see review by Beuthe ^0]). We describe the electron by a wave 
packet, the initial nuclei A and B - by stationary wave functions, while the outgoing particles 
- by plane waves as representatives of a complete set of orthogonal states. 

We assume that the initial nuclei A and B are bound in crystals and their wave functions 
are energy eigenstates localized near central point x^ and xb, respectively, with the uncer- 
tainty of the order of crystal spacing cla,b- Their wave functions are taken as a product: 

^The equality signs in equations throughout the paper should be taken "with a grain of salt" because we omit 
some obvious factors, such as coupling constant, (2ir)~ 3 etc. This makes the formulas easier to read without 
influencing the physical results related to oscillations. 



(2) 
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V a (x) = ^a(x-x a ) e- itE " , (4) 

where a = A,B and E a is the energy of nucleus in the potential of a crystal cell, which is 
equal to the difference of M a and binding energy. The Fourier transform of ip a (x) , which we 
need in what follows, is 

K a (<la) = j *# (x - xje-^™) (5) 

By assumption, nuclei A and B are at rest and thus K(q a ) is centered near q a — with 
uncertainty a a ~ a^ 1 . 

The wave function of the incident electron is taken as a wave packet: 

* e (x) = J dq e K e (q e - p e ) e ^(x-x e )-^ e t ; (6) 

where E e (q e ) = y/q£ + ml , the Fourier amplitude K e (q e — p e ) is centered near q e = p e with 
uncertainty a e , and the maximum of the envelope of the packet is at the point x e at the 
moment t = 0. 

According to the measurements of Novosibirsk group (Pinaev et al., |ll).|12j). which 
used an undulator at electron storage ring, a e > 0.7 x 10 -6 eV = (30 cm) -1 . The theoretical 
analysis of this data and the mechanism of reduction of the wave packet of a relativistic 
charged particle by emission of a photon has been performed by Faleev |13| . (For the general 
theory of wave packets see lectures by Glauber [13].) As for the upper bound on a e , it is 
smaller then W~ 3 E e for an electron accelerator according to the PDG ( |15j . pp. 239-241). 

For the wave functions of the outgoing particles we can take any complete and orthogonal 
set of functions, the most convenient would be just plane waves: 

= e ^=-*Pc* ( 7 ) 

where c = C,D, I. 

It is clear that fermionic nature of neutrino (as well as of e and I) is not essential in the 
problem at high enough energies, therefore we replace the neutrino Green function by the 
Green function of a scalar particle of mass mj, where j enumerates neutrino mass eigenstates, 
j = 1,2,3. Thus: 

Gj(x,t) = - — -r / -= — 5 dkduj . (8) 

n ' (2vr) 4 j w 2 -k 2 - mj 2 + i£ y ' 

We integrate in (j2J) over (t\ + i 2 ) and (i 2 — t±) using 



J 4tt|xi-x 2 | 



and obtain 



Aj = - J dxidx 2 dq e K e (q e - p e )e iqe(xi Xe) V>A(xi ~ x j4 )V'b(x 2 - x B ) 

p ikj\xi— X2| 

g-ipoxx-ipxjxa-iKxa <5(£ e (q e ) + £a - " w) , (10) 

47T Xl — X 2 
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where 



to = Ei + Ed — Eb 



(11) 



and kj = \Ju 2 — m?-. The delta- function in ea. (|lU|) corresponding to the energy conservation 

for the whole process appears due to integration over (t\ + £2)- 

We would like to emphasize that the energy ui of the virtual neutrino does not depend 
on j. The energy uj is determined by energy conservation at point B. Due to stationarity of 
nuclei B,D and lepton I (see eqs.@ and (|7]l) the energy uj is the same for all neutrino mass 
eigenstates. 

Note that e lfc jl x i- x 2 j s the only factor in ea.()lU|) which depends on j and leads to mass- 
dependent effects such as oscillations. Since the initial nuclei A and B are in the potential 
wells, they are well localized near the points xa and x#, respectively; therefore the integrand 
in eq. ljlOj) is essentially different from zero if |xi — x^| < and |x2 — x B \ < as- Hence the 
effective range of integration over xi and X2 is limited by ||xi — X2I — |x^ — x B \\ < a, where 
o-A ~ a-B ~ a- Furthermore, kj differs slightly from uj, so we may expand 



kj = uj + ... 

3 2uj 

In what follows we introduce the following notation: 



(12) 



XAB 



XA - X B 



(13) 



Then 



Jk j \x 1 -x 2 \ _ e ikjX AB exp 




e ikjXAB e ia)(|xi-x 2 |-XAs) 



x 2 



XAB) 



/ ml 

1 + 



(14) 



Lattice spacing a is of order of 10 -8 cm. As for the uj (which is roughly equal to the 
electron energy) , it can vary depending on the beam energy of the accelerator which produces 
electrons. Let us take, for example, uj = 1 GeV. This energy is sufficient for muons to be 
created in the detector. If we take 3 eV as the upper bound for neutrino mass we get 
amj/(2uj) < 10~ 12 , hence the j-dependent correction in brackets in ea.p4jl can be neglected, 
and we obtain from eqs. l|10|) and (|14j) : 



dxidx 2 dc[ e K e (q e - p e )e 



iq e (xi— x e ) e — ip c xi-ip D X2— ipjx 2 



e iuj(\x.i-X2\-XAB) 

i/ja(*i - xa)i/j B (x 2 - x B )—— —5(E e (q e ) + E A - E c - uj) , (15) 

47T|Xi — X2 I 

We see that e lkjXAB is the only j-dependent factor in ea, ()15|) . which means that the phase 
difference between Ai and Aj has the following form : 



4>j = {ki - kj)xAB- 



(16) 



There is no term, containing energy difference, in the right-hand side of eq. (|16|) . This 
results from equality of energies of different neutrino mass eigenstates, emphasized earlier. 
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The form of the phase in ea. (|16|) allows to resolve the problem "equal energy vs equal 
momentum" in favor of equal energy scenario. 

From eqs. ()12|) and (fTK|) we obtain the standard expression: 

<Pij = — 2^~ XAB > ( 17 ) 

where 5m 2 - = mf — m 2 -. 

To integrate over xi and X2 in eq,(|15|) we take into account that a a, as <C xab, and 
hence we can use the expansion 

[Xi -x 2 | ~ x A B ~ n(xi - x A ) + n(x 2 - x B ) , (18) 

where 

n = . (19) 

xab 

From eas. (|18(l and ((To^) we obtain (up to a constant phase factor): 

J-kjXAB r 

A 3 = / dq e e"^ x ^ Xe ^ e (q e - p^Ka^K 3(^)6 {E in - E fin ) (20) 

XAB J 

where 

Ein = Ea + Eb + E e , 
Efin = Ec + Ed + Ei , 
= PC + - q e , 
qs = Pi + Pd - . 

As follows from ea. (|2Uj) and the expressions for the wave functions presented above, the 
amplitude is a function of the vector n, neutrino mass, rrij, energies of the initial nuclei, Ea 
and Eb, central momentum of the incoming electron p e , and momenta of the final particles 

PC, PD, Pi 

In the next section we will use an explicit formula for the amplitude Aj fea. (j2U|0 in the 
case of certain simple expressions for the crystal potential and for the wave packet of the 
electron. We assume that the electron is described by one-dimensional Gaussian wave packet 
with definite direction e = p e /|p e |, that is 

e?q e = 5 ( - e ) \q e \ 2 d\q e \ dQ e (21) 



KeiSie - Pe) = exp 



Pe) 21 



2al 



(22) 



We also assume that the potential is an oscillator near the center of the crystal cell and 
a constant at large distances from this point. Furthermore we consider the initial nuclei in 
the ground states in their crystal cells: 

/Q(q ) = exp (-■&), (23) 



2^' 

where a = A,B. After integration over q e in ea. (|2Uj) by using delta-function in the integrand 
of eq,(|20[) and delta-function in eq.(|21|) we obtain 
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exp 



(q° e e - p e y 
2ol 



where 
and 



Aj = E° e q° e 



Q iqO(x A -x e ) 



exp 



47TXAB 

(am + pc - q°e) 
2a\ 



exp 



(Pi + PD ~ ^m) 
2a% 



E e = Ei + Ec + Ed - Ea — Eb, 



£ = V(£°) 2 -m§. 



Thus introduced coincides with p e = [p e | within cr e . 



(24) 

(25) 
(26) 



3 Integration over phase space 

In this section we integrate the probability of the process (JQ) over the phase volume of the final 
nuclei and over the energy of the final lepton, assuming explicit form of the crystal potential 
and of the electron wave packet (see the end of the previous section) . We demonstrate that 
due to the neutrino energy dispersion the oscillations are suppressed at large distances. We 
think that this result is valid for arbitrary shapes of wave packets. 

We consider a situation when only the direction of the final lepton is measured, while its 
energy is not measured and final nuclei are not registered 2 . Thus we are interested in the 
differential probability for the final lepton to be detected in the solid angle <if2; : 



dP, 



CM, 



f 2 dpc dp D pfdpi _ ^ 



(27) 



Here pi = |p;|; pc, Pd and pi are the momenta of the final nuclei and the final lepton, 
correspondingly, and 



ij 



Pi 



dp c dpp pfdpi 
2Eq 2Ed 2Ei 



where 



Pij — A i Aj . 



(28) 
(29) 



We interpret the electron detection as v e survival, and (jl or r detection as or v T 
appearance, respectively. 
From eq.(|24|) we get: 



exp l t^f-XAB 



J AB 



exp 



(<?°e - Pe 



[un + pc - g° e ) 



(Pi + PD ~ wn) 



a 



B 



■ (30) 



Here we delete inessential factors l/(4-7r) 2 and (E®q®) 2 . In what follows to simplify formulas 
we take 



2 In fact, as it is clear from calculations presented in the Appendix, if we do measure the energy of the final 
lepton, but the precision of the measurement is worse then cr e , the result does not change essentially. 
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M c = M B =M ~ 100 GeV, a A = °b = a ~ 1 keV; 10 MeV < E e < 1 GeV. 
For such choice of parameters the electron is relativistic, while 



^«1, 
M 

where p e = |p e |. 

We calculate Vij, defined by ea.l|28|l. in the Appendix for two cases: 
1. Large a e : 



(31) 



a e » a 



Pe_ 



2. Small a P : 



(32) 



(33) 



We assume that vectors e and n have a substantially nonzero angle between them and the 
module of their difference is of order of unity. We present the results in the form 3 



T> • X AB 

Vij ~ exp It— I exp 



XAB 
j^sup 



where L° s - C = 2E u /5mf^ and, in our case, E v = p e (l + O (ff))- 4 For the first case (|32 



r sup 



(see ea. (|A.30|l ). while for the second case (|3*3*)l 



L 



sup 



roscM 



(34) 



(35) 



(36) 



(see eq.dHHJ). 

The origin of the above suppression is neutrino energy dispersion, on which the oscillation 
length depends. A rather lengthy way to derive eas, ()34j) - ()36|) is given in the Appendix. 

A simple qualitative estimate of the suppression, based on consideration of two particle 
reaction e + A — > v + C, will be presented in ref. |16j. 

It is interesting to note that the suppression length for large a e (the first case) is equal to 
that which arises due to spatial separation of the neutrino wave packets with different mj, 
and hence different velocities, in the case when neutrinos are described not by a propagator, 
but by a wave function (see Dolgov et al. [IB] . Nussinov Kayser fTHi and Dolgov pH]). 

As for the case of vanishingly small a e , there is no separation of neutrino wave packets 
considered in refs.|16j-|19j. but suppression exists due to virtual neutrino energy dispersion. 

Let us end up this section by considering a toy model of two neutrino flavours {y e and 
z^, for definiteness). In this case from eqs. (|27j) and (|34[) we easily find 



dftf, 



1 



-(sin20) 2 



1 



(sin26») 5 



1 — exp 



exp 



X AB 
T2 

■'-'sup 



cos ( XAB \ 
\Lqsc J 



W(\), 



■'sup / 



xab\ 



L 



W(\) , 



(37) 



VI ' 



U e \ = cos 9, U e 2 = sin0 and the factor W(l) does not depend on the lepton 



where 1 

flavour I in the limit of ultra-relativistic final leptons. 



3 The exponential character of suppression arises from the Gaussian form of the initial wave functions in p- 
4 Quite often L OBC is defined as AirE^/Sm 2 . 
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4 Nuclei bound in crystals at non-zero 
temperature 

Now we show that for a non-zero temperature the results do not essentially change. In this 
case we cannot assume nuclei wave functions in crystal cells to be stationary. A general 
expression for the wave functions of the initial nuclei is the following: 

¥A,B(x,t) = ^Cab^bM e- iE Xs\ (38) 

n 

where n enumerates the states with definite energies E\ B , while C\ B is the amplitude of 
probability to measure energy E\ B in the state with wave function ^^(x, t). Using the 
linearity of amplitudes with respect to initial wave functions we obtain 

Pij = ^ ] C m C m iA^ m 'Aj m , (39) 
m rn' 

where for convenience we define the two-dimensional index m = (n\,n2), C m = C"£C r B , 
and Ai m stands for A% calculated with ij/^ and i^g 2 - For typographical reasons there is no 
difference between upper and lower indices. 

Since we consider the case of a thermal equilibrium, we have to average Pij using C m C^, = 
fimm'e~ Em ^ T '> where E m = E 7 ^ + E B 2 . After averaging we obtain 

Pij = y ] e T Ai m Aj m . (40) 

m 

We see that taking temperature into account results in a simple averaging of Pij over 
different stationary initial states weighted with e~ Em ^ T , which is a common case in statistical 
mechanics. Such a correction does not result in any observable effect due to the smallness of 
thermal energy T compared with uj. Indeed, 

5m% ( 5m% 4 T \ 

<t>ij = -^xab + O —^-x AB , 41 

2,UJ \ 2,UJ UJ J 

where uj corresponds to the case of zero temperature (ground states, tt-2 = 0): 

uj = E 1 + E d -E%. (42) 

We see that the correction is essential if xab ~ (u/T)L osc ~ 10 10 L osc for u = 1 GeV, 
T = 300 K. In the previous section it was shown that at such distances oscillations are washed 
out. Thus we may safely neglect this correction and recover to the standard expression (|17|). 



5 Conclusions 

l.Our calculations explicitly confirm the equal energy scenario, advocated by Lipkin and 
Stodolsky (see text after eqs, ((TTj) and (|16jl). 

2. The oscillation phase has its standard form, see ea.(|17|). 

3. When integrated over the phase volume of the final nuclei and over the energy of the final 
lepton, the oscillating term in the probability of the process under consideration vanishes 
exponentially at the distances L sup = (2p e /a e )L osc for large a e , and L sup ~ (M/a)L°? c for 
small a e including plane wave limit for electron; see eqs. (|33|) and (|3f)|) . 

4. For non-vanishing temperature T the standard expression for the phase difference is valid 
up to correction ~ (T/uj). 
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7 Appendix A 

1. General formulas 

The aim of this Appendix is to integrate Pij given by eq. (|30j) over the phase space of 
nuclei C and D and over the energy of the final lepton /: 



v* = I r,*£*£«^. (A.i) 

3 1 3 2E C 2E D 2 



In what follows for simplicity we consider an ultra-relativistic electron: 

q°e = E° e , (A.2) 

where E® and are defined by eqs, (|25j) and (|2T)|) , We assume that the target and the 
detector contain the same nuclei: 

M A = M B = M c = M D = M, a A = a B = a. (A.3) 

The Appendix consists of six parts. In part 1 feas. (|A.l|) - (|A.13|0 we make no special 
assumptions about the width <j e of the electron wave packet. In parts 2 and 3 (eqs. (|A.14|) - 
()A.30|) ) the case of a relatively broad electron wave packet is considered: 

a e » (A.4) 
The parts 4 and 5 (eqs, (|A,31|) - (jA.41|0 are devoted to the narrow electron wave packet: 

a e « (A.5) 
Part 6 (eqs. (|A.42|) - (|A.44|l ) deals with the case of an electron plane wave: 

a e = 0. (A.6) 
Taking into account eq. (|A.3|) and omitting pre-exponential factor in eq. (|3()j) . we get: 

5vn?;XAB \ 



Pij = exp i- 



2lu 



(E° e e - Pe ) 2 (am + Pc - E°e) 2 (con - P/ - p D ] 



exp \ — v\ ^ ** (A '* ! 

where E® is defined by ea. (|23|) . 
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As uj does not depend 5 on pc, it is convenient to integrate Pij in ea. ()A.l|) first over pc 
and then over p£, and E[. For this purpose we introduce function I\\ 



I\ e l exp 



(um -pi- pp) 2 \ dp D pidEi 
a 2 J 2E D 2 



where 



h 



exp 



(£ e u e - Pe ) 2 (am + p c - E° e e 



0^2 



Here 6 is the angle between pc and -E^e — un, while 

dpc 



pc sin 6* (f# dEc ■ 



2E C 



7T sinOdOpcdEc- 



(A. 



(A.9) 



(A.10) 



In this Appendix pc and denote modules of three- vectors pc and pr>, respectively. 
Let us integrate over 9 first: 



/ 



dE c 



\E° e 



e — wn 



■ exp 



(E° e - Pt 



\2l 



err 



exp 



\E°e 



Lun\ 



PC? 



cxp 



(|^°e-wn| 



(A.ll) 



It is evident that l\ is not singular at zero angle a(e, n) between neutrino and electron. 
Still we assume that vectors e and n have a substantially nonzero angle between them and 
the module of their difference is of order of unity, that is 



E% — wnl » cr, 



(A.12) 



and thus we may neglect the second term in the curly brackets in eq. ([A.ll|) . 

As was mentioned in the body of the text, the suppression of oscillations is a consequence 
of dispersion of neutrino energy. For a(e, n) = and large a e the suppression length is given 
by ea. (|3oT) because electron energy spread coincides with that of neutrino. 

If we consider now the case of very small a e , the suppression governed by nuclear a 
becomes dominating. Let us assume, for simplicity, that m e = and nucleus C does not 
interact with the crystal. For zero angle production of neutrino the momentum of the nucleus 
C is equal to zero up to a, while its energy does not exceed the value of the order of a 2 /M, 



5 We could have defined ui not by ea. 1)11(1 for vertex B, but by energy conservation at vertex A. However, the 
oscillation phase still would not depend on momentum of nucleus C . 



w{q° e ,p c ) = E c (pc) ~ E e (q° e ) - E A , 



<7g depends on pc'- 



in such a way that 



7e=<Zc°(Pc) 



M£(Pc),Pc) dE c {pc) dE e (q° e ) 



dp, 



c 



dp c 



dq® dp, 



= 0. 



c 
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and hence neutrino energy equals E e up to this value. We omit the derivation noting only 
that the suppression length in this case is very large: 



j sup rose PeM 
ij ~ V 1 ' 



a- 



(A.13) 



2. Large cr e , integration over pc- 



Now we consider the case of relatively broad electron wave packet. We assume that the 
first exponent in the curly brackets is much sharper than one which is out of the brackets, 
and cuts out the effective region of integration. The condition under which one can make 
such assumptions will be obtained further (see eq. (|A.17|) ). 

The dominant exponent (the first in the curly brackets) has a maximum when 



f(E c ) = (\E° e (E c )e - un\ - pc{E c )f = 0. 



(A.14) 



Note that E® is "chosen" by global conservation of energy (ea.(|25j)). and hence depends on 
Ec- 

Let us denote the solution of this equation as E c . It has the sense of the most probable 
value of E c . Evidently p° c = J {E c ) 2 - M 2 . For non- vanishing values of the angle a(e, n) 



defined above p c ~ |p e — um| ~ p e . Here and in what follows the sign "~" means "is of the 
order of magnitude" . 



We expand f(Ec) up to the third order in (Ec — E c ): 

f(Ec) = \{E C - Ecf^riEc) + \{Ec- &C?^V%) 
Taking into account that p c ~ p e we estimate the derivatives: 



(A.15) 



1 d 2 f(E c ) 

2 dE 2 c 



M_, 

v n - 

Pc 



p\ 



1 d 3 f(E c ) M 3 



6 dEl 



{P C Y 



M 3 



(A.16) 



Now we can see that the effective width of the exponent in the curly brackets in (jA.ll|) is 
ap e /M while the width of the first exponent in (|A.llj) is a e . Thus the electron wave packet 
may be considered broad one if 



Pe_ 

M' 



(A.17) 



From eqs. ()A.ll|) . (|A.15|) and ()A.16|) we have 

^M 3 {E C -E C ) 3 {E C -E C ) 2 



h ^ / dE c 



exp 



l^e - urn 



exp 



p c ra 2 



or 



(E c + u;-E A -p e ) 2 _ fap c Y 



exp 



\a e Mj 

M 2 (e c -E c + (E c + u- Ea — p e ) (^) 2 ) 



{ap c f 



(A.18) 
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Note that the integrand in I\ is not vanishingly small if 



\E c -E A -p e 



< 



(A.19) 



E C -E° c + (E° c + u; - E A - Pe ) 



~ M 



(A.20) 



and thus taking into account eq. (|A.17|) we may estimate the essential range of Ec — E c : 



\Ec-E c \<—. 



(A.21) 



This means that the exponent in the first line of ea. (IA.19|) may be omitted. Furthermore, 

the small corrections in the second and third lines of ea, ()A,19|) proportional to are 
negligible. Using this we obtain 



h 



dE c - 



|£°e - un 
and, after integration: 



exp 



(7 



(E° c + u;-E A - p e f M\E C - E° c ) 2 



/l = M eXP 



(E° c -p e -E A + uj) 



(A.22) 



(A.23) 



3. Large a e ; integration over ~pD and E\. 

First we integrate eq. (|A.8|) over the angle between pd and urn — pi and obtain: 



Pi 



dE D dEi- 



Pi 



11 ^ ,,xp 



|wn - p/| 
(|wn - p;| -pd) 21 



. ° m ij x AB 

-e 1 2^ exp 



(E° c -p e -E A + uj) 2 



of. 



exp 



(\un - pi \ +pn) 



(7 



Assuming that the angle between vectors n and 1 = pi/pi satisfies the inequality 



a 



a(n,l) > — , 

Ve 



we disregard the second exponent in the curly brackets in ea. (|A.24|) . 
To integrate over Ed and E\ we introduce 

F(E D , Ei) = X 1 {E c -p e -E A + uj) 2 + (|um - p,l| - p D f- 



(A.24) 



(A.25) 



(A.26) 



Expanding all functions of Ed and E\ in the exponents in eq. (|A.24|) around the most 
probable values E D and Ef, which are solution of the equation 



F(E D ,E l ) = 0, 



(A.27) 



we obtain from eq, (|A,24|) 



13 



Vi 



Pi 



o 5 



- E u n\ M 



exp 



M 2 



(AE D ) 2 - 2( 



1 



2E„ 



AAf 



dEpdEie 



-(AE D +AEi) 



-)&E D t\E l -(— + —)(m) 



or 



(A.28) 



where AE D = E D - E° D , AE t = E t - Ef, 



(^n-p°l)(n-^-l) 
\E„n - pPlj ' 



(A.29) 



Pi=\J (Ei) 2 — mf, mi stands for the mass of the final lepton, p® D = </ (E^) 2 — M 2 and 
^ = E° D + Ef - E B . 

Taking into account that A ~ 1 we may integrate over Ejj and E\ in eq. (|A.28j) : 



<r> %PeOe 
Pij = o -775- exp 



M 2 



2E„ 



IE,, IE,, 



where we use eqs. (fA~26jl and (fA~27j) which define E° D and E<j>. 
4. Small er e ; integration over pc- 



(A.30) 



Let us now consider the case of small o e , that is 

Pe 

The integration will be mainly carried out as for the case of large o e . 
The first exponent in ea. (|A.ll|) has a sharp maximum when 

E c = E° c =p e + E A + E B -Ei- E D . 



(A.31) 



(A.32) 



As before, Eq has the sense of the most probable value of Ec- We would like to emphasize 
that the definition of Eq here is different from one which is for the case of large o e (see 
eas. (jA.14|) and (|A.32|l ). To estimate the contribution of the first exponent in the curly 
brackets in eq. (|A.ll|) it is convenient to use the function /(Eq) defined in eq. (|A.14|) . 

Expanding f(Ec) near E c to the first order of AEq = Eq — E c we obtain for I\ from 
ea. (|A~TT|) : 



\p e e - un\ 



dE c exp( ^ — ) exp 



or 



'■>C 



Since 



— (E r ) ~ a — 
dE c K C) Pe 



then the term containing derivative in ea. (|A.33|) is negligible, 

—z—-{E c )t\h c r < 1, 

o 2 dE c o p e 



(A.33) 



(A.34) 



(A.35) 
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and can be omitted. Granting this we integrate in ()A.33|) and obtain 

(jgeg - un\ -p° c ) 2 ' 



a 2 a e 



\p e e - tun 



exp 



(A.36) 



where p^ = y {E^) 2 — M 2 . Note that exponents in eas. (|A.23|) and (|A.36|) are different. The 
value of a e determines which exponent in ea. ()A.ll|) is more narrow and cuts out the effective 
region of integration. The wider exponent remains after integration in (|A.11|) . 



5. Small a e ; integration over po and E\. 

Similarly to the case of large a e integration in 
tun — pi gives: 



over the angle between pjj and 



dEidE D - 



Pi<y^v e 



exp 



13 J " \Pe£ — om| • I pi — am I 

Sm^XAB (\p e e-um\ - p° c ) 2 (|pj-am| - p D ) 2 



2uj 



(7 



(A.37) 



where we take into account the assumption (|A.25|) . 

We integrate over Ed and E[ analogously to the case of large a e . Introducing 



G(E D ,Ei) = (\p e e - om| - p° c ) 2 + (\p t - um\ - p D 



(A.38) 



and expanding all functions of Ed and E[ in exponents in ea. (|A.37|) near the most probable 
values E D and Ef, which are defined by equation 



G(E D ,E l )=0, 



(A.39) 



we obtain 



V, 



4 
Pl<T(T e 
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exp 



\p e e - E u i\\ ■ \pfl - E v x\ 



■ e 2E V 



dE D dEie 2 ( E -) 



5m ■ ■ x a r 



M 2 



(7+ 



(AE D ) 2 (AE D ) 2 2AE D AE l (AE t ) 



+ 



+ 



+ 



(p D ) 2 p 2 (e — n) 2 p 2 {e — n) 2 p 2 (e — n) 2 



(A.40) 



where AE D = E D - E° D , AEi = E x - E?, p? = y ' {Ef) 2 - mf, p° D = ^(E° D ) 2 - M 2 and 
E V = E° D + Ef-E B . 

After integration in (|A.4fl|) taking into account eqs. ()A.38|) and (|A.39|) we finally obtain 



Vij a ~M 2 6XP 



dmjjXAB ( Smf-XAB p e \e - n\a" 



2E„ 



2Eu 



2MEL 



(A.41) 



6. The case of electron plane wave. 



In this part of Appendix we would like to discuss separately the case of the incident 
electron with definite momentum p e . As is well known in the limil a e — > Gaussian 
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(2703/2 0-3 ex P ( (qe CT f e) ) -> <5 3 (q e - p e ). By substituting it in eq.® #e(q e - Pe) =5(q e - Pe 
we immediately obtain the following expression for the amplitude 



ikjXAB JVe (xa — x«) 



(«1a) 



(A.42) 



where = Pc + wn ~~ Pe 5 q_B = Pi + Pd — <^ n ; an d as previously we assume the nuclei to 
be in the oscillator-like potential. 
We use 

(6{E in - Efin)) 2 = T 5(E in - E fin ), (A.43) 

where T is a duration of the experiment, to square the delta-function in eq. (|A.42|) . and obtain 
the probability normalized per unit of time: 



Pij = a 2 5(E in - Ef in ) exp 



. Smjj (um + Pc ~ Pe) 2 (pi+Pd- wn) 
i— - — x AB k ^ 



2u 



o- 



a 



(A.44) 



We may integrate eq. (|A.44|l over pc, Pd and E\ exactly as in the case of small a e (see 
eqs. (|A37jt - fPnjn . This proced ure results in ea. ()A.41|) . 
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